/ 9 9^0/ '7V7 2 _ 


N94 


Scattering by Cavity-backed Antennas on a 

Circular Cylinder 

Leo C. Kempel and John L. Volakis 
Radiation Laboratory 
1301 Beal Ave. 

Ann Arbor, MI 48109-2122 

November 15, 1993 


j ?&/ ¥/ 

35 

21888 


Abstract 

Conformal arrays are popular antennas for aircraft, spacecraft and 
land vehicle platforms due to their inherent low weight and drag prop- 
erties. However, to date there has been a dearth of rigorous analyt- 
ical or numerical solutions to aid the designer. In fact, it has been 
common practice to use limited measurements and planar approxi- 
mations in designing such non-planar antennas. In this paper, we 
extend the finite element-boundary integral method to scattering by 
cavity-backed structures in an infinite, metallic cylinder. In particular, 
we discuss the formulation specifics such as weight functions, dyadic 
Green’s function, implementation details and particular difficulties in- 
herent to cylindrical structures. Special care is taken to ensure that 
the resulting computer program has low memory demand and mini- 
mal computational requirements. Scattering results are presented and 
validated as much as possible. 
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1 Introduction 

Conformal antenna arrays are attractive for aircraft, spacecraft, and land 
vehicle applications since these antenna systems have low weight, low drag, 
flexibility, and cost advantages over conventional protruding antennas. The 
majority of previous studies pertaining to non-planar conformal antennas 
has been conducted experimentally due to a dearth of rigorous analysis tech- 
niques. Traditional rigorous techniques involve an integral equation and are 
limited in terms of radius of curvature and structural complexity. Some ap- 
proximate methods have been considered but these are restricted in accuracy 
and element shape. 

Recently, the finite element-boundary integral (FEM-BI) method was 
successfully employed for the analysis of large cavity-backed planar arrays 
[1], The resulting system is sparse due to the local nature of the finite ele- 
ment method, whereas the boundary integral sub-matrix is fully populated. 
However, by resorting to an iterative solver such as the Biconjugate Gradi- 
ent (BiCG) method, the boundary integral sub-system may be cast in circu- 
lant form allowing use of the Fast Fourier Transform (FFT) in performing 
the matrix-vector products. This BiCG-FFT solution scheme ensures 0(N) 
memory demand for the entire FEM-BI system and minimizes the computa- 
tional requirements. 

In this paper, the FEM-BI formulation is extended to scattering by aper- 
ture antennas conformal to a cylindrical metallic surface. In contrast to 
the planar aperture array, the implementation of the cylindrically conformal 
array requires shell-shaped elements rather than bricks, and the required ex- 
ternal Green’s function must satisfy the boundary conditions on J.he surface 
of the cylinder. In its exact form, this Green’s function is an infinite series 
which imposes unacceptable computational burdens on the method. How- 
ever, for large radius cylinders, a suitable asymptotic formula is available 
and herein used for an efficient evaluation oTthe Green’s function. In ad- 
dition, the resulting BI system is again cast in circulant form to ensure an 
0(N) memory demand and to take advantage of the FFT’s efficiency when 
carrying out the matrix-vector product. 

A primary difficulty in studying cavity-backed antennas mounted on curved 
surfaces is the lack of reference data. In this paper, scattering calculations 
based on the FEM-BI method are compared with data based on different 
techniques. Although such validation is necessarily limited, it provides con- 
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fidence in the formulation’s accuracy so that this approach may be used in 
extending the available reference data. 


2 FEM-BI for Circular Cylinders 

In this section, the FEM-BI formulation is developed for cavities recessed 
in an infinite metallic cylinder, having walls which coincide with constant 
p-,4>- or z-surfaces (see figure 1). As usual, the finite element formulation 
permits substantial modeling flexibility, including cavity inhomogeneities, 
lumped loads, super/substrate antenna configurations, or microstrip line and 
so on. 

The FEM-BI approach possesses both low memory and computational 
demand when implemented with a BiCG-FFT solver. Although the system 
of equations associated with the FEM formulation is sparse, the boundary 
integral sub-matrix is fully populated. However, if the aperture mesh is a 
uniform grid, the BiCG-FFT solver may be employed for that portion of 
the system thus retaining (9(N) memory demand for the entire system. In 
addition, the solver require low computational demand since both sparse 
matrix- vector products and discrete convolutions using FFTs require only 
O(NlogN) operations per iteration. 

The FEM-BI formulation begins with the weak form of the vector wave 
equation followed by specification of appropriate vector shape functions and 
dyadic Green’s function. The resulting FEM-BI equations are then used to 
solve for the total electric fields within the cavity and on the aperture (see 
for example Volakis et al. [2]). For the specific configuration at hand, the 
weak form of the wave equation can be written as 

r f v x • v x Wj(p, z) 

Jv, \ p r {p,<f>,z) 

-kle r (p,<f>,z)Wj(p,</),z) • Wi(p, <f>,z)^pdpd<f)dz 

+{k 0 a) 2 8 a (j)6 a (i) I I [Wi(a,<f>,z)- p(a,<f),z)x 

J 5 | J Sj 

V 2 {a,j>,z) x p(a,<f> , z) ■ Wj(a, (j> , z')] dj dz d<f> dz = f\ nt + f' xt (1) 
In this, Wi are vector basis functions with support over the volume Vi which is 
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Figure 1: Illustration of a typical cavity-backed antenna situated on a metal- 
lic cylinder and the associated coordinate system. 
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associated with the i th degree of freedom, and in a similar fashion, S, and 
represent aperture surfaces associated with the i th and j th degrees of freedom, 
respectively. The appropriate dyadic Green’s function is denoted by G 2 and it 
has convolutional = — <f>,z = z — z) form when evaluated on the surface 

of the cylinder, p = a. The unprimed coordinates represent the test point 
while the primed ones denote the source point. The free-space propagation 
constant is given by fc 0 = where A 0 is the free-space wavelength. The 
cavity is filled with an inhomogeneous material having relative constitutive 
properties e r and p T . The function S a (t)S a (j) is the product of two Kronecker 
delta functions. Hence, it identifies which pairs of unknowns belong to the 
aperture and accordingly contribute to the boundary integral sub-matrix. 

The FEM-BI equation (1) may be rewritten in matrix form as 



[ 0 ] [ 0 ] 
[ 0 ] [ 0 ] 


[*?} _ un 

~ { 0 } 


(2) 


where the entries of [*4] are due to the FEM portion of the formulation and 
[&] is the boundary integral sub-matrix. In (2), Ej P and £j nt denote degrees 
of freedom associated with the aperture and interior fields, respectively. In 
this, ff xt are functions of the external excitation and will be discussed later 
in the paper. 

An important factor in choosing the finite elements for gridding the cav- 
ity is the element’s suitability for satisfying the mathematical requirements 
of the formulation as well as the physical features of the antenna system. 
Traditional node-based finite elements associate the degrees of freedom with 
the nodal fields and have proven unsatisfactory for three-dimensional electro- 
magnetics applications since they do not correctly represent the null space 
of the curl operator and hence spurious modes are generated [3, 4]. In con- 
trast, edge- based elements correctly model the curl operator and therefore 
the electromagnetic fields. In addition, edge-based elements avoid explicit 
specification of the fields at corners where edge conditions may require a sin- 
gularity. Jin and Volakis [6] presented edge-based brick elements which are 
convenient for rectangular-type structures and cavities. For cavities residing 
in a circular cylinder, shell elements are the natural choice. 

Cylindrical shell elements possess both geometrical fidelity and simplicity 
for cylindrical-rectangular cavities. Figure 2 illustrates a typical shell element 
which has eight nodes connected by twelve edges: four edges aligned along 
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each of the three orthogonal directions of the cylindrical coordinate system. 
Each element is associated with twelve vector shape functions given by 


Wn (p, 4>i z) = W p (p, 4>, 2 ; •, 4 > r , z t , +), 
^*h&iPi 4*i z) = Hp(p, 4*i Z, '1 4*r 1 Zb, )i 

H^14 {Pi 4*i z) H /7 <^ ( P 1 4*i Z, Pbi '1 Zfi *4" )? 
Wss (p, 4*1 z) = WfiPi 4>, 2 ; p b , ■, z b , -), 


W&{p,4>,z) = W p (p.4>,z-,-,4>i,z t , -) 
Hy(p, 4*1 z) = Hyp, 0 , 2 ; •, 4> h z b , +) 

H^sCPi 4*i z) H / ^(p, 2, p a , •, 2j, ) 

W 67 (p, 4*1 Z^ H^^p, 4*,Z , Pa ? '1 ^fci “I”) 


'W'x s(/?i ^) “ -2TJ Pbi y 4“ )? W^26 (Pj *0 Pen ^ r > *J ) 

h / 4s(p, 4*1 z) = Hyp, 4>i z-, pb, 4>i, ■ 1 -), H / 37 (p, 4 , 2 ) = Hyp, <f>, 2 ; p a , 4>i , •, +) 

(3) 

where Hy is associated with the edge which is delimited by local nodes (l,k) 
as shown in figure 2. As seen from (3), three fundamental vector weight 
functions are required for the complete representation of the shell element. 
They are 


W p (p,4*,z;p, 4*i z,s) 
H ’^{p,4*,z;p,4>,z,s) 
W 2 (p,4>,z;p, 4*i z,s) 


S£b (4 - <ft)(z - z ) . 

ah p 

^{p - p)(z - z)4* 

y-ip- p){4*~ 4*)z 


(4) 


where the element parameters (p a , pi>, <f>^ (f> r , z\> y z t ) are shown in figure 2, t = 
Pb — p a , a = 4> r — <f>i and h = z t — z^ The ^-term which appears in the 
definition of the ^-directed weight (4) is essential in satisfying the divergence 
free requirement, i.e. so that V • Wj — 0 1 . Note that as the radius of the 
cylinder becomes large, the curvature of these elements decreases, resulting 
in weight functions which are functionally similar to the bricks presented 

1 Wj(p,<f>,z) will only satisfy this requirement within the volume of the element. These 
weighting functions introduce artificial charges on the faces of the element and are not di- 
vergenceless at element interfaces. This is allowable since these elements do not guarantee 
normal field continuity across the element faces. 
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by Jin and Volakis [6]. Having specified the vector basis functions, we may 
proceed to develop the matrix entries for the system (2). 

The FEM-BI system is composed of two parts: a sparse FEM matrix and 
a fully populated BI sub-matrix as shown in (2). The FEM matrix entries 
are represented by 


A - — 

— *st 


Hr 


- k l f r 

K o t r 1 st 


(2)0 


(5) 


where constant material properties have been assumed within each element. 
The subscripts (i,j) refer to the row and column of the matrix entry and 
correspond to the test and source edges, respectively. The auxiliary functions 


V x W t {p,<}>,z\p x ,<t> x ,z t ,Si)pdpd<}>dz 

1 st ~ jy Pjs Zjs <j) x z\ pi, <j> xx Zi, S{)p dp d<f> dz (6) 

are identically zero unless both test and source edges share at least one el- 
ement in common, resulting in a highly sparse system. Physically, such a 
system is a consequence of the locality property inherent in a partial differ- 
ential equation formulation. In (6), the direction of the source and test edges 
are represented by ( s,t ) <E {/>, respectively. Since the edges of the mesh 
are aligned along three orthogonal directions, only six combinations of ( s,t ) 
are required for and only three such combinations for and all of these 
are evaluated in Appendix A. Since (6) is symmetric with respect to source 
and test edges, the FEM matrix will also be symmetric. 

A lumped impedance post may be included in the formulation by adding a 
term to (1) and equivalently to (5); surface or sub-surface metallization layers 
may also be modeled. Radially oriented lumped loads are approximated 
in the FEM-BI formulation by a filamentary load located at ( 4>l,zl ) [2]. 
Such posts have length /, cross-sectional area s and impedance Zl . The 
contribution to [.4] is given by 


Aij = jkoZo-1- f - ^ ZL) Wi{p, <f>, z)Wj(p x <f>, z)pdpd<f>dz(7) 

SZi Jv p 

which may be readily evaluated in closed form. In addition, infinitesimally 
thin metallization layers may be represented by simply fixing a priori the 
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weight coefficients to zero for weights associated with edges which are tan- 
gential to the metal. This is a consequence of using a total electric field 
formulation. The symmetry and sparsity of the FEM system [«4] is main- 
tained after the addition of these loads while the BI system \Q) remains fully 
populated and symmetric. 

The boundary integral provides an exact boundary condition for mesh 
closure and its construction relies on a cylindrical dyadic Green’s function. 
The entries of the boundary integral sub-matrix are 


Gij = {k 0 a) 2 / / W t {a, <f>, z\ pi, <j>i, £,, s t ) ■ 
Js, Js, 

p{a, <f>, z ) x (? 2 (a, 4>, 2 ) x p{a, <f> , *') 

• M a ( u <, (f) z , pj , (f)j , Zj , s j ) d<f> d d(f) dz 


( 8 ) 


where the weight functions are given by (4) and evaluated at the surface 
p = a. In (8), the dyadic Green’s function ((*2) satisfies both the radiation 
condition and the Neumann boundary condition at p = a. This dyadic 
Green’s function may be expressed exactly [8] 


G**(a,<j>,z) — - J7r-yi E I" (¥) ~ ~£t 7) e j(n *~ kti) dk z 

_L.fr ^ } (7) e Hn*- k ,* )dk 

^ 2) (7) 


G*’(a,lz) = 


r*+(a j _ 1 V r I ^ 2) (7) (nk z V H ™ I 

( 27r ) 2 7 [tfi 2) (7) U<> afc P/ H' n {2) 


H (2) i (7) 

(7) 


e j(n4>-krz) dk 


(9) 


where 7 = A^a and k p = y — A;^. However, for large radius cylinders, (9) 
is computationally prohibitive. In these cases, which are of main concern 
in this paper, it is advantageous to employ an asymptotic expression for G 2 
[9, 10, 11, 12]. These employ a creeping wave series expansion of which only 
the two direct path contributions (see figure 3) are retained. The formula 
due to Pathak and Wang [9] 

G z *(a,4>,z) ~ -^qe~ 3k ° a ^(cos 2 0 + q(l - q)(2 - 3cos 2 9)'j v(/3)^ 
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C?^*(a, <f>, c) ~ }k °" sin 0 cos (1 — 3g(l — q)) n(/?) | 

G*(a, 4>, z) ~ ( 57n ^ + 9(1 - ?)( 2 - 3sm 2 0)) t>(/?) 

+9 [ 5 ec 2 ^(u(/9)- y (/3))]| 


(10) 


" 2 A i • 

where fi = ks ^, -^j 3 and q = -^ has proven quite accurate. In the defi- 
nition of /?, s is the usual geodesic path length (s = y (a<f>) 2 + z 2 ) and 6 is 
the direction of the trajectory = tan~ l ). Depending on which of the 

two direct paths (shown in figure 3) is used, ^ or $ = 27r — <j>. The soft 
and hard Fock functions, u(fi) and v{j3) respectively, are characteristic of 
on-surface creeping wave interactions and have been extensively investigated 
by Logan [13]. Although computation of the Green’s function (10) is now 
tractable, evaluation of (8) must be done so that a discrete convolutional 
system is maintained and the singularity of (8) at s = 0 is properly treated. 

Care must be taken in evaluating (8) so that the overall storage require- 
ment remains 0(N) and the singular integrals of (8) are accurately computed. 
If uniform zoning is used, the resulting sub-matrix ([£/]) is block Toeplitz and 
hence amenable to solution using the BiCG-FFT method. For the non-selfcell 
contributions, mid-point integration may be used while a regularization pro- 
cedure must be employed for the self-cell. Bird [12] noted that (10) recovers 
the metallic screen Green’s function when /? = 0 within the available ap- 
proximation order. This suggests that (8) may be regularized by adding and 
subtracting from (10) the function 


2(j 0 ( fl , = 


= VV 

I + ~V~ 

K o J 


g J koR 

27 7r 


R=\f-r\ 


(11) 


which is the free-space dyadic Green’s function multiplied by two. The re- 
sulting regularized Green’s function (curvature contribution) is given by 


G 2 


! (a,<f>,z) ~ — jfco3 | (cos 2 6 + <jr(l — q)(2 — 3cos 2 9f) [u(/?) — 1] 
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G U {a,ftz) 


jk 0 

2x 


qe jkoS sin 9 cos (1 — 3^(1 — q)) [t>(/?) — 1] j 
— j (sm 2 0 + q( 1 — q)( 2 — 3s77? 2 0)j [u(/3) — 1] 

«( 0 ) -»(/*))] } 


+q $ec 2 0(i 


(12) 


and since it is no longer singular it may be evaluated numerically. The planar 
contribution may be calculated in the manner described previously by Jin 
and Volakis [5]. The FEM-BI matrix has now been fully developed and it 
remains to specify the excitation function for external sources. 


3 Plane Wave Excitation 

Plane wave excitation of the geometry is considered in this section for scatter- 
ing analysis. The use of the exact boundary condition in (1) allows coupling 
of an exterior excitation field into the cavity. We will describe the form of 
the source functional, / t ext , and discuss its numerical implementation. 

The forcing functional, due to exterior sources (/, ?xt ) is given by 

fi Xt — }Z 0 k 0 a [ <f>' , z) • /5(a, ft , z) x H^(a,ft y z')dft dz' (13) 

Js, 

where W{(p, <f>, z) is the testing weight for the i th row of the matrix and H cyl 
represents the magnetic field on the cylinder’s surface in the absence of the 
cavity. A plane wave 

E { = 

W = Y 0 (ic ‘ x 

= Y 0 [p* sin 7 cos Oi - ft cos 7 -z sin 7 sin 0,] e jk ° [psine ‘ cos( *-* l)+z cosfl d 

(14) 

is assumed to be incident on the cylinder from the direction (0,-,^,) where 
7 is the polarization angle and e' = 0'cos 7 + ^‘sin 7 is the electric field 
polarization. The total surface field is given by the sum of the incident 
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and corresponding scattered field from the infinite metallic cylinder [14]. 
Specifically, 


+ (a, *,=) 

= iHf + zHf 


(15) 


where 


= -2H 


0 jk 0 cos 6 t z oo 


cos 7 


nk 0 a sin 0, n f^ 00 [//< 2) (*; 0 fl sin 0.) 
j n sin 7 cos 1 cin( - +0 _ 0i) 

/: 0 a sin 6, Hn 2 \k 0 a sin#,). 


H?{a^z) = 


i2r 0 ^-e jfcoCO8fl<2 T 

nk 0 a 




Ti— “oo [Hn^\k 0 a sin #,) J 


(16) 


is obtained from traditional modal analysis. These expressions may be ap- 
proximated by retaining only a few terms of the series if k 0 as\n8i is small. 
However, as this parameter becomes large (e.g. for large a and #, — ► 90°), 
(16) may be replaced with equivalent asymptotic representations similar to 
those considered earlier. Utilizing Watson’s transformation and Fock theory 
[14] in connection with (16), we find that 


Tjcyl 

*1$ ~ 


Hf ~ —Y 0 sin 7 sin 0 ie i k ° coa<> < z e- jk ° a8 ' m6 '*’’ [<7 (0) (m$ p )]* 

P= 1 

2 2 

j2Y 0 cos 7 - ™ e ikoCoa8 ' z £ [/<°>( m $ p )] * 

0 * p=l 

2 r 

—Y 0 sin 7 cos #,e jtoC08 *' 2 £(-l y> e -H**W, g M(m$ p ) 

p=i L 

(17) 


in which j - {<f> - <f>i), $ 2 = (4> ~ <f>i ) ~ f , m = and 

complex conjugation is denoted by an asterisk. The appropriate iar-zone 
Fock functions (g^°\ g^ and f^) are given by Logan [13]. 
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The asymptotic formulas (17) are quite accurate except in the geometrical 
optics region (<J> % <^,). In this case, Goriainov’s [15] expressions 


H? 


have been found to be more accurate and can be used instead of (17). The 
Fock functions (G, G ^ and F ) are again defined in Logan [13]. These 
surface field expressions may be used to calculate the entries of the column 
vector { ff xt } efficiently via a numerical evaluation of (13). In particular, 
the modal series (16) is used when & o asin0, < 10 and either (17) or (18) 
for k 0 as\n6i > 10 as appropriate. With the excitation functional and the 
FEM-BI matrix now specified, the BiCG-FFT method [16, 17] may be used 
to determine the unknown electric fields within the cavity. 

4 Scattering 

Once the cavity aperture and volume electric fields have been determined by 
solving (2) for an external excitation, the radar cross section (RCS) may be 
calculated. The far-zone fields may be computed by integrating the aperture 
fields with a suitable Green’s function. In this section we present the relevant 


—Y 0 sin o sin 0 i e , * oC “ e * 2 je -, *° asin e ‘* 1 ^ (0) (m$,)]* 

+eJ fc 0 a.infl, co. (*-*,) | G (_ mcos (^ _ | 

j2Y 0 cosa - e i*oCc^z| e -ifc.a.inft», [/(O^j)]* 

+e jk 0 a .in 0 , co. (*-*) [F(-mCOs(^- &))]* J 

+ y; sin a cos 0i e jk ° cos * j e ~ jk ° a sin e ' *• 5 (0) ( m 4> , ) 

b; 9(,,(m * i) ] 

^(7(_ m cos (^ — <^ t )) 


— e 


—J 


m 

k 0 a sin I 

jk Q a sin 8% cos (4>—4 


m 


k 0 a sin 0 X 


G<»( — m cos ( <(> 


-«n 


(18) 
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formula for calculating the far-zone fields and hence the RCS due to excitation 
by a plane wave (14). 

To determine the far-zone fields, we begin with the integral representation 
for the scattered magnetic field in terms of the aperture fields. We have 

H’(r,6,<t>) = jY 0 k 0 a f V 2 (r,0,<t>]a,(f>' ,z') - 

J s 

jp(a, <f > , z ) x E(a, 4 > , z )] d<t> dz (19) 

with (r,0,<f>) indicating the observation point in spherical coordinates. When 
the observation point is very far from the cylinder, the dyadic Green’s func- 
tion in (19) can be replaced by its far-zone representation 

—jk Q r 

Z7a(r, 6, 4>- a , <j>\ z) ~ C — - \G h 0$ -f G Bz 0z -fi G+*W] (20) 

where the unprimed unit vectors are functions of the observation position 
and the primed ones are functions of the integration point in (19). The 
components of this far-zone Green’s function 


G h ~ 

r^Bz 

& ~ 
G M ~ 


l 2 kp COS 0 jkoCOS g z ' ^ n 

(2tt) 2 ( k Q a sin 0) 2 H' n {2 \k 0 a sin 0) 

•9-00 i , 

3 1 c jk 0 cos Bz ST' e jn(^+(0-<£ )) 

(2-k) 2 a n=-oo Hn 2 \k o asin0) 

9 - OQ 1 , 

i i JkoCosOt y' i e in(Z +(4>-4 )) 

(27r) 2 a sin 0 H (2 \k Q a sin 0) 


( 21 ) 


are determined by a mode matching procedure. As one might expect, these 
series converge rather slowly for large k Q a sin 0. They must therefore be recast 
in another form by employing Watson’s transformation and Fock theory as 
was done previously (17). In doing so, we obtain 


G h 


G 0Z ~ 


G++ ~ 


koJX>S0 e jk 0 cos 9z Y^-iy e ~ ik ^(°)(m$ p ) - i-^^5 (1) ( m$ p) 


k o sm0 0 j J- 




j k 0 cos Bz g —j k 0 a sin 9$ p 

p=l 


m 1 * * 

1 ~ jkpCotiBz ST 


2air sin 0 


p=I 




( 22 ) 
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where the Fock functions are the same as those used with (17) due to reci- 
procity. As was the case for the plane wave source, Goriainov’s [15] approx- 
imations are more accurate in the geometrical optics region (<p s; <f>) and 
similar expressions may be obtained for (22) as was found for (17). The far- 
zone scattered field can be computed numerically by using (19) and either 
the series or asymptotic formula as appropriate. Having done so, the RCS is 
calculated from 




lim 47rr 2 

r— ► oo 


\H‘( r,M)l 


(23) 


Above we presented a FEM-BI formulation suitable for modeling cavity- 
backed structures embedded in a circular cylinder. Next, we consider a few 
numerical calculations aimed at validating this formulation and in giving 
us an appreciation on how the cylinder’s curvature influences the scattering 
parameters. 


5 Results 

Having solved for the electric fields induced by an incident plane wave, the 
resulting RCS data must be validated with known results. As previously 
mentioned, available measured or computed data is rather scarce and as 
a consequence, we are forced to rely on limiting cases in order to validate 
this work. As the radius of curvature decreases, a cylindrical-rectangular 
cavity will approximate a planar-rectangular cavity. Another limiting case 
involves comparison of an elongated 3-D cavity with a corresponding 2-D 
cavity for normal incidence (0, = 90°). Finally, we may compare our infinite 
cylinder results with a finite Body of Revolution (BOR) model for certain 
polarizations and angles of incidence. We begin with the quasi-planar case. 

The first validation effort for scattering by cavity-backed patch antennas 
relies on the fact that a small patch on a very large radius cylinder is quasi- 
planar and approximates rather well an equal sized planar patch. For our test 
we chose as a reference a planar 1.448” x 1.083” patch residing on a 2.89" x 
2.10 x 0.057" cavity filled with a dielectric having e r = 4. The equivalent 
patch on a 10A cylinder is 6.46° x 1.083" residing on a 12.90° x 2. 1 0” x 0.057" 
cavity. Figure 4 shows the results for the patch on a large radius cylinder 
with corresponding data for the planar cavity-backed patch. Clearly, the two 
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Monostatic RCS (o/X ) [dB] 



Angle (0) [deg] 


Figure 4: Comparison of RCS for a planar patch (1.488" x 1.083") residing on 
a 2.89 x 2.10 x 0.059 cavity filled with e r = 4 dielectric and a corresponding 
quasi-planar patch on a large radius (IOAq) cylinder. 
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RCS patterns are in excellent agreement, and although figure 4 illustrates 
only monostatic scattering in the <f> = 0° plane, additional runs for normally 
incident monostatic scattering and various bistatic situations yield similar 
agreement. 

Comparisons may also be made for elongated cavities and 2-D MoM re- 
sults. Long narrow cavities have very little axial interaction for principal 
plane ( 6 = 90°) excitation and therefore results based on this formulation 
should compare well with corresponding 2-D data. It is well known that the 
RCS of the 3-D scattering body of length L Ao is related to the corre- 
sponding 2-D scattering of the same cross section via the relation 



Such a comparison is shown in figure 5 for monostatic scattering by a 45° x 
5A x 0.1 A cavity for both principal polarizations. Once again the agreement 
between the two results is excellent, thus providing a partial validation of the 
formulation for highly curved geometries. We remark that similar agreement 
has been observed for bistatic scattering in the 6 — 90° plane. 

The planar approximation eliminates the effects of curvature, which is a 
primary interest in this work, and the 2-D comparisons done above are only 
valid for normal incidence. To consider oblique incidence on a highly curved 
structure, we resort to comparisons with a Body of Revolution (BOR) code 
for wraparound cavities. Since the BOR code can only model finite struc- 
tures, we simulate an infinite cylinder by coherently subtracting the far- zone 
fields of the finite structure without a cavity from similar data which includes 
the cavity. Such an procedure mimics common measurement practices and 
was found suitable for near normal incidence and quite acceptable near graz- 
ing incidence in the case of H-polarization (a = 90°). An example calculation 
for the latter case is given in figure 6 where a bistatic scattering pattern is pre- 
sented in the <f> = 0° plane due to a plane wave incident at (0, = 90°, <^>, = 0°). 
Clearly, there is good agreement between the FEM-BI results and data based 
on the BOR formulation. 

The previous comparisons serve to validate the formulation. Having done 
so, it is instructive to examine the effect that curvature has on the scattering 
properties of cavity-backed patch antennas. Consider a 2 cm x 3 cm patch 
residing on a 5.0 cm x 6.0 cm x 0.07874 cm cavity which is filled with 
a dielectric having e r = 2.17. Figures 7 and 8 illustrate the behavior of 
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Monostatic RCS (o/A, 2 ) [dB] 



Angle (<|>) [deg] 


Figure 5: Comparison of 2-D MoM results and FEM-BI RCS results for a 
45° x 5A 0 x 0.1 Aq air- filled cavity. 
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Bistatic RCS (a/Xf) [dB] 



Observation Angle (0) [deg] 


Figure 6: Comparison of the RCS computed via the FEM-BI method and a 
BOR code for a 3Ao xO.lAo air-filled wraparound cavity excited by a normally 
incident H-polarized (a = 90°) plane wave. 
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this geometry as a function of frequency and curvature. Evidently, the 
resonance behavior of this patch is sensitive to curvature for both principal 
polarizations. The frequency response for E-polarization is more sensitive to 
curvature since the radiating surface field component is parallel to the long 
side of the patch and cavity. If the patch and cavity were oriented so that the 
long side is in the ^-direction, the response to H-polarization would exhibit 
greater sensitivity. Such an effect is important to low observable antenna 
designers since they want to operate the antenna in the region of lowest 
RCS. This low return region is a consequence of delicate cancellations due to 
the physical layout of the aperture. Such cancellations are not as complete 
for highly curved structures as they are for planar cavities. 

Conformal antenna designers often use wraparound antennas to achieve 
omnidirectional coverage. Two different configurations are typically used: a 
wraparound cavity where the cavity is filled with a single continuous collar 
of dielectric and discrete cavities symmetrically placed around the circum- 
ference of the cylinder. These two configurations are shown in figure (9). 
Since near resonance, the radiation properties of these two types of antennas 
is identical, any RCS advantage which one might possess could govern the 
appropriate choice of arrays. Figure 10 compares the E-polarized monos- 
tatic scattering at 3 GHz in the 6 = 90° plane for a wraparound cavity and 
four discrete cavities; where the patches and cavities are identical to those 
used in the previous example. Not surprisingly, the wraparound structure 
has a higher return due to coupling within the substrate. However, since in 
this case the scattered field is due to the z component of the surface field 
(^-directed magnetic currents), both cavities yield large scattered fields in 
the four directional lobes. Figure 11 is the corresponding comparison for 
H-polarization. In this case, the scattered field is attributed to the <j> compo- 
nent of the surface fields (z-directed magnetic currents). Therefore, substrate 
modes diffract near the patch resulting in discrete lobes for the discrete array 
while creeping waves shed isotropically for the continuous wraparound cav- 
ity. Low observable designs will favor discrete cavity arrays over wraparound 
cavities since the scattering may be channeled in preferred directions and 
the overall scattering level is consistently lower. A final example is shown in 
figure 12 where we observe that other than the expected higher scattering 
from the wraparound cavity, the scattering behavior of the two arrays is very 
similar. 
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RCS (c/T) [dB] 


Frequency [GHz] 


Figure 7: RCS frequency response for a 2 cm x 3 cm x patch residing in a 
5 cm x 6 cm x 0.07874 cm cavity with e r = 2.17 as a function of curvature 
for E-polarization (a = 0°). 




RCS (oA ) [dB] 


Frequency [GHz] 


Figure 8: RCS frequency response for a 2 cm x 3 cm x patch residing in a 
5 cm x 6 cm x 0.07874 cm cavity with e T = 2.17 as a function of curvature 
for H-polarization (a = 90°). 





Figure 9: Illustration of two types of arrays: (a) wraparound array; (b) 
discrete cavity array 
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Monostatic RCS [dBsm] 



Angle (<j>) [deg] 


Figure 10: Comparison of E-polarized monostatic RCS at 3 GHz for a four 
patch array placed on a wraparound collar or in four discrete cavities. The 
patches and cavities are identical to the one used in figure 7. The observation 
plane is 6 = 90°. 
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Monostatic RCS [dBsm] 



Angle (<J>) [deg] 


Figure 11: Comparison of H-polarized monostatic RCS at 3 GHz for a four 
patch array placed on a wraparound collar or in four discrete cavities. The 
patches and cavities are identical to the one used in figure 8. The observation 
plane is 6 = 90°. 
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Monostatic RCS [dBsm] 



Angle (0) [deg] 


Figure 12: Comparison of H-polarized monostatic scattering at 3 GHz by a 
four patch array placed on a wraparound collar or in four discrete cavities. 
The patches and cavities are identical to the one used in figure 8. The 
observation plane is <f> = 0°. 
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6 Conclusions 


In this paper, we have presented a Finite Element Method - Boundary Inte- 
gral technique suitable for electromagnetic scattering calculations for cavities 
embedded within a circular, metallic cylinder. This formulation is analogous 
to the FEM-BI approach used by Jin and Volakis [1, 5, 6] and may accord- 
ingly be used for the analysis of scattering by a large array of cavity-backed 
patch antennas. These cavities need not be identical, periodically spaced or 
homogeneously filled and may in fact may possess lumped impedance loads 
or surface metallization layers. The FEM approach employs vector finite el- 
ements which properly represent the electromagnetic fields and possess high 
geometrical fidelity for cylindrical-rectangular cavities. Such elements were 
presented and are analogous to the bricks used for modeling rectangular 
cavities. In addition, we presented an efficient method for evaluating the on- 
surface and far-zone dyadic Green’s functions. The presented formulation is 
amenable to solution using the BiCG-FFT method provided uniform zoning 
is used across the aperture and as a consequence, this implementation has 
low computational and memory demand. We have presented some validation 
of this work with appropriate limiting cases which provides further archival 
reference data. In addition, we showed how this formulation may be used to 
influence conformal antenna designs. 
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Appendices 


A FEM Matrix Entries 

The matrix entries for the FEM portion of the system (5) are given in this 
appendix assuming that the cylindrical shell elements (4) are used in (6). 
These integrals are given by 

I ( pz = — ^ —J (<f> - <i>s){<f> - <f>t)d<f> 

ot j 

= jiKy (i - p °) + 3 + P 1 ) (A ~Pb) + 2^ aPi _ rtO) + 

(pb - pi) - 2< {ps + pt) + PsPt In J z ( 2 - 2 *)( 2 - 

J* = “7? J p b ( p - p»)(p - P<) d P 

= J^y [ a (2 “ 1 (p* + ^*) + P‘P * ln ) + 

\{pI~pI) J* r (<j>-4>,)(<t>-4> t )dt] 

[ (z - z 3 )(z - z t )dz 

4? = [j (pi -/>;)+! if . + m y. - 4 ) + \pa yt - pi) ] x 
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f {<f> - - <Pt)d<p ( A - 1 ) 

J<Pi 

Each of the above unevaluated integrals is of the form 

j" (£ - £.) ({ - l) di='- (L 2 - V 2 ) ({, + £,) + 5 (t/ 3 - L 3 ) + £,£, (U - L) 

The integrals are used in the assembly of the FEM portion ([.4]) of 

the system. 


B Fock Functions 

The asymptotic form of the dyadic Green’s function with observation both on 
the surface of the cylinder and in the far field involves Fock functions. These 
have been extensively studied and tabulated by Logan [13]. The numerical 
evaluation of these functions are performed either for small arguments or 
large arguments. 

The on-surface Fock functions used in this paper are 

u(0 = f°° — tj— |e -J ^ T dr 

2 V 7T ./ooe-jW 3 W 2 (t) 


,(£) = f 

y/TT Jooe ~> 2lr / 3 W2(t) 


(B-l) 


where w 2 (t) and its derivative w' 2 (t ) denote Airy functions of the Second 
Kind. For small arguments (£ < 0.6), the asymptotic expansion of (B-l) is 
given by 


»(£) ~ 1.0-^fi+j^( 3 + 5^v^e-5£l + ... 

ti(f) ~ 1.0- + + 

while a rapidly converging residue series is used for £ > 0.6 


(B-2) 


10 , 

>(0 ~ = e_jiT " 
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u(0 ~ =2e^^J2( T nV'e~ j * Tn (B-3) 

n=l 


where r„ and r n are zeros of u>2( r ) and w' 2 (t ), respectively. Those zeros are 
given in the following table 


Table B-l 

Zeros of the u> 2 (t) and iu 2 (t) 
T n = |r„|c" J ? and r' n = 

n 

Trt 

\<\ 

i 

2.33811 

1.011879 

2 

4.08795 

3.24819 

3 

5.52056 

4.82010 

4 

6.78661 

6.16331 

5 

7.94413 

7.37218 

6 

9.02265 

8.48849 

7 

10.0402 

9.53545 

8 

11.0085 

10.5277 

9 

11.9300 

11.4751 

10 

12.8288 

12.3848 


The far-zone Fock functions are given by 


* w (0 

/»(« 

G (,) (0 
F (l \ 0 


j l f e* T 

y/n Jr u^(r) 

jl r e it-r 

\fn Jr wi(t ) 

9 0) (O^ 


dr 

dr 


(B-4) 


where w\{t ) and its derivative inj(r) denote Airy functions of the First Kind 
and the integration contour is given by Logan [13]. These functions, 
and /^(£), may be calculated using 
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9 {0) (0 = 


2.0e' j £ £ < —1.3 


c(m) 


= 1.39937 + V — 1-3 < 4 < 0.5 

■ ^ m I 


m=l 


m: 


10 


= 


Kc/frn)^ 


r', c/(m).Ai(m) 
1.8325e 


0.5 < £ < 4.0 


-(0.8823-j0.5094)«-iV] ^ > 4 Q 


(B-5) 


<J (1 ’(f) = £<-2. 


.0^25 


? ) 


2.8 


= ±s^p {ir -^ - 2.8 <e <0.5 

m = l 


10 e [«a(m)^] 


0.5 <{<4.0 


= k y: 4 ./ \ 

T^i ^*( m ) 

= -1.8325 (0.8823 - j0.5094 + j£ 2 ) J-* 0 - 8823 -' 0 - 5 ™ 4 *"^] ^ > 4 0 

(B— 6) 


/ (0) (£) = j% 


. n , /, 0.25 0.5 \ _,fi , 

— + — Je J 3 £<-1.1 


6 c(m) 


= 0.77582 + e~ jir/3 V] -±-f- «) m - 1.1 < £ < 0.5 


= e 


„T, m! 

m=l 

10 

“•” /3 E-n77-r 0.5 <<<4.0 

m = 


= 0.0 £>4.0 


(B-7) 


with constant k = e -J W 6 and the coefficients for (B-5) and (B-6) given in 
the following table 
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Table B-2 

Constants for (B-5) and (B-6) 

m 

c(m) 

a’(m) 

;4z(m) 


0.7473831 

1.01879297 

0.5356566 


-0.6862081 

3.2481975 

-0.41901548 


-2.9495325 

4.82009921 

0.38040647 


-3.4827075 

6.16330736 

-0.35790794 


8.9378967 

7.37217726 

0.34230124 

6 

56.1946214 

8.48848673 

-0.33047623 



9.53544905 

0.32102229 



10.52766040 

-0.31318539 



11.47505663 

0.30651729 



12.38478837 

-0.30073083 


The corresponding constants for (B-7) are given as 


Table B-3 

Constants for (B-7) 

m 

c(m) 


KSBSM 

i 

1.146730417 

2.33810741 

0.70121082 

2 

0.86284558 

4.08794944 

-0.80311137 

3 

-2.0192636 

5.52055983 

0.86520403 

4 

-9.977776 

6.78670809 

-0.91085074 

5 

-14.59904 

7.94413359 

0.94733571 

6 

49.0751 

9.02265085 

-0.97792281 



10.04017434 

1.00437012 



11.00852430 

-1.02773869 



11.93601556 

1.04872065 

10 


12.82877675 

-1.06779386 
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